BRAF mutation predicts survival after immunotherapy across multiple 


cancer types 


Author list: 
Weiting Ge’, Wen Cai'*, Dehao Wu", Wangxiong Hu’, Weidong Han?, Shu 
Zheng*', Hanguang Hu*14 


Weitinge Ge, Wen Cai, and Dehao Wu contributed equally to this work. 


1 Cancer Institute (Key Laboratory of Cancer Prevention and Intervention, 
China National Ministry of Education), the Second Affiliated Hospital, School 
of Medicine, Zhejiang University, Hangzhou, China 

2 Department of Gastroenterology, the Second Affiliated Hospital, School of 
Medicine, Zhejiang University, Hangzhou, China 

3 Department of Medical Oncology, Biomedical Research Center, Sir Run Run 
Shaw Hospital, School of Medicine, Zhejiang University, Hangzhou, China 

4 Department of Medical Oncology, the Second Affiliated Hospital, School of 


Medicine, Zhejiang University, Hangzhou, China 


Full address for correspondence: 

Hanguang Hu; Department of Medical Oncology the Second Affiliated 
Hospital, School of Medicine, Zhejiang University, Hangzhou, China; E-mail: 
huhanguang@zju.edu.cn; Tel: 86-571- 87784718 


1 


Shu Zheng; Cancer Institute (Key Laboratory of Cancer Prevention and 
Intervention, China National Ministry of Education), the Second Affiliated 
Hospital, School of Medicine, Zhejiang University, 88 Jiefang Road, 


Hangzhou, China; E-mail: zhengshu@zju.edu.cn; Tel: 86-571-87783768 


Abstract: 

Recently studies in selected tumors suggested that BRAF mutation may 
associates with survival benefit from immune checkpoint inhibitor (ICI) therapy" 
3. To broadly investigate this association at a pan-cancer level, we analyzed two 
independent ICI treatment cohorts (MSKCC*: n = 1630, and Dana-Farber>: n = 
249). BRAF-mutant patients exhibit better overall survival in the MSKCC cohort 
(Hazard ratio [HR] = 0.55, 95% confidence interval [Cl]: 0.43-0.72; P<.001) and 
the result is validated by the Dana-Farber cohort (HR = 0.68, 95% Cl: 0.46-0.99; 
P = .045). A multivariate analysis adjusting tumor mutational burden, mismatch 
repair status, cancer type, age and sex confirmed the results (adjusted HR = 
0.58, 95% Cl = 0.43-0.78; P < .001). Immunogenomic features analysis of 
TCGA dataset indicated that patients may respond to immunotherapy in various 
mechanisms. This finding substantially improve the therapeutic prospects for a 


sizeable fraction of patients who benefit from immunotherapy. 


Immune-checkpoint inhibitor (ICI) therapy has revolutionized the 
treatment of advanced cancer. However, durable benefit is limited to a 
minority of patients. Previously studies have suggested that PD-L1 
immunochemistry, deficient mismatch repair (AMMR), tumor mutational 
burden (TMB) may correlate with clinical response*®. BRAF mutation 
associates with bad prognosis and chemotherapy resistance in most of 
cancers, which new therapy choices were urgently needed. Recently, a phase 
Ill trial reported that BRAF-mutant melanoma patients who received first-line 
ICI treatment showed better overall survivals than wild type patients?. Also, 
BRAF-mutant glioblastoma and non-small cell lung cancer patients reported 
potential benefits from ICI treatment in retrospective studies'*. However, it is 
unclear whether BRAF mutation associates with ICI treatment outcome 
across diverse human cancers. In this pan-cancer study, we identified that 
BRAF-mutant tumors associate with a better outcome after ICI treatment and 
present a small number of tumor neoantigens. Our study also indicates that 
patients with/without BRAF mutation may respond to immunotherapy in 


various mechanisms. 


Methods 
All clinical and mutational data including TMB score were obtained from 
the cBioPortal database (https://www.cioportal.org) on May 25, 2020. The 


accession ID of two ICI treatment cohorts are TMB_MSKCC_ 20184 and 


MIXED_ALLEN_2018°. The immunogenomic data and MMR status were 
obtained from the published analyses of in the cohorts of TCGA PanCancer 
Atlas’ and MSKCC-IMPACT>, respectively. All cases without BRAF mutation 
status or cancer type information were excluded. The Kaplan—Meier survival 
analysis with a two-sided log-rank test was employed to estimate BRAF 
mutation in relation to overall survival (OS). A multivariate Cox regression 
analysis was performed to determine the contribution of mutation to OS. A 
two-sided Mann-Whitney U test was used to examine the differences between 
groups. Statistical analysis was conducted with R v3.6.1. All analyses were 
conducted using deidentified public-use data. This study was deemed exempt 
from institutional review board approval and patient informed consent was 


waived. 


Results 

After excluded 8518 cases without BRAF mutation status or cancer type 
information, we obtained 38142 patients across 30 cancer types in curated set 
of non-redundant studies from cBioPortal. Totally, 1292 patients with BRAF 
V600E, 858 patients with other BRAF mutations, and an overall mutation 
frequency of 5.6% were identified (Fig.1). The top three cancer types with the 
highest frequency of BRAF mutation are thyroid cancer (53.9%), melanoma 
(45.0%) and colorectal cancer (13.6%). 

We further employed the MSKCC ICI treatment cohort* to investigate the 
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association between BRAF mutation and OS. A Kaplan—Meier survival 
analysis was performed on 1630 patients with follow up data, including ten 
types of advanced cancer as previously described*. As showed in Fig2a, 
BRAF-mutant patients exhibit a significantly better OS than wild type (Hazard 
ratio [HR] = 0.55, 95% confidence interval [Cl]: 0.43-0.72; P < .001). This 
result is validated by an independent ICI treatment cohort of Dana-Farber 
Cancer Institute including 249 patients with seven kinds of cancer (HR = 0.68, 
95% Cl: 0.46-0.99; P = .045; Fig2b). To adjust potential bias, we selected 
1078 cases with sufficient data including BRAF mutation, TMB score, MMR 
status, cancer type, age and sex from the MSKCC cohort and perform a 
multivariate Cox regression analysis. After adjusting the above factors, the 
association between BRAF mutation and improved overall survival after ICI 
treatment is statistically significant (adjusted HR = 0.58, 95% Cl = 0.43-0.78; 
P < .001). The patient's clinical characteristics of two cohorts were showed in 
Supplementary Table 1 and 2. 

We further employed TCGA immunogenomic data’ (n = 7251) to 
investigate the potential mechanism which could explain immunotherapy 
efficacy of BRAF-mutant patients. Although an increased number of tumor 
neoantigens has been linked to better response during immunotherapy’, 
patients without BRAF mutation exhibit a significantly higher number of 
neoantigens (P < 0.001, Fig3a) rather than BRAF-mutant group. Then we 
compared intratumor heterogeneity (ITH) between groups. Previous studies 
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stated that increased ITH associates with lower efficacy of immunotherapy in 
multiple cancer types’. As showed in Fig3b, BRAF-mutant patients exhibit a 
significantly lower ITH than the wild type patients (P < .001). We also 
compared six immune subtypes in order to systematically evaluate the 
immune landscape of cancer’. The results showed that the C3 (Inflammatory) 
subtype is enriched in BRAF-mutant group (55.4% versus 21.4%), whereas 
the C4 (lymphocyte depleted) subtype is relatively higher in wild type group 


(13.3% versus 6.1%, Fig3c). 


Discussion 

Our study firstly reported that BRAF-mutant patients who received ICI 
treatment presented a better overall survival than wild type at a pan-cancer 
level. This finding will improve the situation that clinical response of ICI 
treatment is limited to a minority of unselected patients. Higher TMB score 
and dMMR, which have been identified as positive predictive markers*?, 
define about 20% patients of MSKCC ICI treatment cohort. Our findings will 
include BRAF-mutant patients and substantially improve the therapeutic 
prospects for a sizeable fraction of patients (28.7%). We still need to identify 
better biomarkers to define the likelihood of benefit from immunotherapy. 

Higher TMB or dMMR patients occupy a higher number of tumor 
neoantigens presented on MHC molecules. ICI treatment promotes an anti- 
tumor immune response to recognize these neoantigens as foreign. Our study 
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found that BRAF-mutant tumors present a small number of tumor 
neoantigens. However, these patients show a better clinical outcome after ICI 
treatment. This result suggests that mechanisms of immunotherapy response 
in BRAF-mutant tumors may not be entirely determined by tumor neoantigen 
presentation. Lower ITH is one of mechanisms to improve immunotherapy 
efficacy®. Additionally, an immune landscape analysis shows that the C3 
subtype, which is characterized by a balanced immune response’, occupies 
more than a half of BRAF-mutant tumors. In contrast, low lymphocyte 
infiltration subtype (C4) accounted for small fraction in BRAF-mutant group, 
less than half of wild type group. Lower ITH and immune equilibrium that 
featured in our study partly explained good response of BRAF-mutant patients 
after ICI treatment. Therefore the relationship of BRAF and tumor immune 
microenvironment as previously reported"? deserves further analysis. 

There are some limitations of our study. First, the findings of this study 
were mainly based on MSKCC and Dana-Farber ICI treatment cohorts, which 
didn’t contain all kinds of cancer with frequently mutated BRAF (eg: thyroid 
cancer). Therefore, BRAF mutation still has great potential to identify more 
patients who response to immunotherapy in various cancer types. Second, 
our study made a bioinformatics analysis to investigated potential 
mechanisms of BRAF-mutant immunotherapy, which can inspire perspective 
studies and need to be validated in clinical trials. Also, the specific 


mechanisms of immunotherapy response of tumors with different driver 


mutation or various microenvironment are not fully explained and worth to be 


further investigated. 
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Figure legend 

Figure 1. Prevalence of BRAF mutation in 38142 patients across 30 cancer 
types. GIST: gastrointestinal stromal tumor; CNS: central nervous system. 
Figure 2. Kaplan—Meier survival analysis of BRAF mutation and overall 
survival after immunotherapy. Patients were obtained from MSKCC ICI 
treatment cohort and Dana-Farber ICI treatment cohort. ICI: immune 
checkpoint inhibitor; HR: hazard ratio; Mut: mutation; WT: wild type. 

Figure 3. Immunogenomic differences of patients with/without BRAF mutation. 


Patients were obtained from TCGA PanCancer Altas. 
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Figure 2. 
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Figure 3. 
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Supplementary Table 1. Patients’ clinical characteristics in Mixed_Allen_2018 cohort 


Wild type BRAF mut P value 
(n=168) (n=81) 
Age*, mean +SD 61.8 +13.0 56.1 +17.0 0.023 
Sex, n (%) 0.760 
Fmale 63 (37.50) 32 (39.51) 
Male 105 (62.50) 49 (60.49) 
Cancer Type, n (%) <0.001 


Melanoma 74 (44.05) 77 (95.06) 
Non-Small Cell Lung Cancer 53 (31.55) 3 (3.70 
26 ( 


(3.70) 
Bladder Cancer 15.48) 1 (1.23) 
Head and Neck Cancer 12 (7.14) 0 (0.00) 
Anal Cancer 1 (0.60) 0 (0.00) 
Small Cell Lung Cancer 1 (0.60) 0 (0.00) 
Soft Tissue Sarcoma 1 (0.60) 0 (0.00) 

Drug Type, n (%) <0.001 
anti-CTLA-4 72 (42.86) 73 (90.12) 
anti-PD-1/PD-L1 90 (53.57) 4 (4.94) 

Combo 6 (3.57) 4 (4.94) 


*Data on age was missing for 66 patients of the wild type group, and 31 patients of the 
BRAF-mutant group. Statistical analyses were performed using the Statistical Package 
for the Social Science (SPSS), version 20. Difference in age was analyzed with t-test. 
Difference in sex was analyzed with Chi-square test. Difference in cancer type and 


drug type were analyzed with Fisher exact test. 
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Supplementary Table 2. Patients’ clinical characteristics in TMB MSKCC 2018 cohort 


Wild type BRAF mut P value 
(n=1461) (n=169) 
Age, mean +SD 61.6+13.4 60.5+15.3 0.319 
0.821 
Fmale 549 (37.58) 62 (36.69) 
Male 912 (62.42) 107 (63.31) 

Cancer Type, n (%) <0.001 

Melanoma 199 (13.62) 118 (69.82) 
Non-Small Cell Lung Cancer 311 (21.29) 24 (14.20) 
Colorectal Cancer 97 (6.64) 13 (7.69) 
Bladder Cancer 206 (14.10) 6 (3.55) 
Cancer of Unknown Primary 82 (5.61) 4 (2.37) 
Head and Neck Cancer 136 (9.31) 2 (1.18) 
Esophagogastric Cancer 122 (8.35) 1 (0.59) 
Glioma 114 (7.80) 1 (0.59) 
Renal Cell Carcinoma 150 (10.27) 0 (0.00) 
Breast Cancer 43 (2.94) 0 (0.00) 

Skin Cancer, Non-Melanoma 1 (0.07) 0 (0.00) 

Drug Type, n (%) <0.001 
anti-CTLA4 70 (4.79) 29 (17.16) 
anti-PD-1/PD-L1 1193 (81.66) 85 (50.30) 

Combo 198 (13.55) 55 (32.54) 


*Statistical analyses were performed using the Statistical Package for the 
Social Science (SPSS), version 20. Difference in age was analyzed with t-test. 
Difference in sex and drug type were analyzed with Chi-square test. Difference 


in cancer type was analyzed with Fisher exact test 
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